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We study the renormalized energy-momentum tensor of cosmological scalar fluctuations during 
the slow-rollover regime for power-law inflation and find that it is characterized by a negative 
energy density at the leading order, with the same time behaviour as the background energy. The 
average expansion rate appears decreased by the back-reaction of the effective energy of cosmological 
fluctuations, but this value is comparable with the energy of background only if inflation starts at 
a Planckian energy. We also find that, for this particular model, the first and second order inflaton 
fluctuations are decoupled and satisfy the same equation of motion. To conclude, the fourth order 
adiabatic expansion for the inflaton scalar field is evaluated for a general potential V(4>)- 
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I. INTRODUCTION 
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£3 , 

Inflationary models (see [l], 0] for a textbook review) , born to solve the problems associated with the standard big- 
bang theory, have their greatest success in the generation of the initial condition for the large scale structure of the 
Universe through the amplification of the quantum cosmological fluctuations during their accelerated era [|| . Using 
the last cosmic microwave background (CMB) data ([J]) and considering the transition from inflation to the standard 
Big Bang cosmology the constraints on the amplitude and spectrum of CMB fluctuation become constraints on the 
qq ■ different models of inflation. All these constraints depend on the linear treatment of cosmological perturbations, 
so, on using the most recent data, we are close to measuring the details of the spectrum of fluctuations. What is 
i-H i still not completely understood is their energy content and the eventual back-reaction on the inflationary expansion 
OA responsible for their amplification. 

Among the interesting effects, appreciable only beyond the linear order, one notes, for example, the effect of the 
non-Gaussianities in the matter power spectrum and in CMB anisotropics (see [5] for a review). However, from the 
theoretical point of view, we feel that the back-reaction of gravitational fluctuations on thegeometry is one of the most 
interesting issues @. Within the inflationary context, this problem has been tackled in II] with the intringuing 
y| result that the energy-momentum tensor (EMT henceforth) of fluctuations may slow down inflation. Thanks to this 
5h ■ result in the last decade there has been a renewed interest in the subject of back-reaction [9lJlCillll.ll2l.ll3l.ll4l.ll5l.ll6|. 

The aim of this paper is to compute the renormalized EMT of cosmological fluctuations for a pow er-law model of 
. inflation, according to the adiabatic regularization scheme [13] already used in previous papers jl3l. [T§|. While the 
adiabatic vacuum for a cosmological scalar fluctuations, and its associated EMT, can be computed for generic space- 
times, the unrenormalized EMT can be calculated analytically only if exact analytic solutions for the field Fourier 
modes are available. This is the case for the power-law inflation [19[ where analytic solutions for a scalar field are 
available. The main theoretical problem which arises is that inflation never comes to an end because of its slow-roll 
parameter e is constant. 

Power-law solutions can be also applied in different backgrounds as, for example, the alternative cosmological 
scenario called Ekpyrotic Universe [2Q|, and in the Pre-Big Bang scenario (for a review see [2l|) in string cosmology. 

The plan of the paper is as follow. In Sec. [TT] we present the linear cosmological perturbations in the uniform 
curvature gauge (UCG henceforth). In Sec. IIIII we extend the UCG to second order and using the Einstein equations 
to second order we give the equation of motion for the second order inflaton fluctuation. In Sec. HVI we give the solutions 
for first order fluctuations for the particular case of power-law inflation. In Sec. |V]we compute the renormalized value 
for the correlator. We discuss the EMT and the back-reaction on the geometry in Sec. I VII and we give our conclusions 
in Sec. IVIII In the Appendix we exhibit the fourth order adiabatic expansion for a general potential V((p). 
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II. LINEAR PERTURBATIONS IN THE UNIFORM CURVATURE GAUGE 



We consider inflation in a flat universe driven by a classical minimally coupled scalar field with a general potential 
V(<p). The action is: 



d i xC 



d A x, 
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where C is the Lagrangian density. 

Let us study the fluctuations of the scalar field ip(t, x) around its homogeneous classical value <fi(t) and include 
metric perturbations. For the homogeneous case we have 



ZH<j> + V> = , 
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(2) 
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where H = a/a is the Hubble parameter and a is the scale factor, using the notation = 1/(8%G) for the (reduced) 
Planck mass definition. 

The scalar perturbations around a spatially flat Robertson- Walker metric are (ds 2 — g^dx^dx"): 
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where Dij = didj — 1/3 V 2 <5y. In the above we neglect vector and tensor perturbations. If we choose to work in the 
UCG one obtains: 



ds 



-(1 + 2a)dt 2 - a/3, i dtdx i + a 2 5 ij dx i dx j 



(5) 



We note that this choice fixes uniquely the gauge, just as the more frequently used longitudinal gauge (for a review of 
cosmological perturbations in this gauge see [22]), The equation of motion for the first order scalar field fluctuations 
in the UCG is given by [131 ] : 
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where a dot denotes a derivative w.r.t. the time t. Using the energy and momentum constraints in their linearized 
version 
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The Fourier tranform modes of v = atp satisfy: 
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where ' denotes a derivative with respect to the conformal time ?y, drj = dt/a. On comparing the last equation with 
Eq. (12) of [23[ it is immediate to see that ip satisfies the same equation as the Mukhanov variable Q. Therefore, the 
uniform curvature gauge has the advantage of singling out the true dynamical degrees of freedom (the matter ones) , 
even if it has the disadvantage of being non diagonal in the metric perturbations. 



III. BEYOND THE LINEAR ORDER 



To second order we consider a metric having the following coefficients: 

500 = -1 - 2a - 2a (2) 
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The above metric is the extension of the uniform curvature gauge to second order: a*- 2 ) and [3^ are the scalar 
perturbations to second order. To second order, scalar, vector and tensor perturbations do not evolve independently 
as is the case in first order. For this reason we take into account second order vector and tensor perturbations, 



(2,') i 

represented by the divergenceless vector x) an d by the transverse and traceless tensor h\*' , respectively. In the 
above we have omitted first order vector perturbations (which die away kinematically) and tensor perturbations 
(which satisfy the usual equation h + 3Hh — V 2 h/a 2 = 0). With this approximation we are neglecting the EMT of 
vector and tensor perturbations, and their correlations with the scalar perturbations. We finally note that the choice 
in Eq. (llip (including vector and tensor metric elements to first order) fixes the gauge completely to second order. 

Using the energy and momentum contrains up to second order we can obtain the equations of motion for , 0^ 
and <p®, in particular one obtains for the second order field fluctuation the following equation (see fi~3l |) 
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where s is a non-local spatial contribution given by 
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IV. POWER-LAW SOLUTION FOR LINEAR PERTURBATION 

We now want restrict ourselves to the case of power- law inflation where a(t) ~ t p with p > 1. In this case we 
have an exponential potential 



V(</>) = V exp 



A 



(14) 



with V = ^P(3p - 1) and A = (§) 



1/2 



The scale factor and the omogeneus solution are given by 
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respectively, and the slow-roll parameters by 
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The equation of motion ([9]) has an exact solution for this particular case [19|. In Fourier space, substituting the 
effective mass squared = p-(3p — 1) and the Hubble factor H — p/t, one obtains 



k 2 

ip k + 3Htp k + —tp k = . 

Introducing ipk — a 3 ^ 2 (fik and the conformal time rj the equation becomes 
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which gives the general solution 
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with v = | + ^ry- We now consider quantized fluctuations of the inflaton: 



<p(r],x.) = 



(2tt) 3 / 2 
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where the b k are time-independent Heisenberg operators. In order to have the usual commutation relations among 
the fej.: 



[6 k , UA = [i>l i>l] = o [k, hi] = s& (k - k') 

one must normalize the solution to the equations of motion through the Wronskian condition: 



(19) 



V>k <f' k - (p' k <P* k = — 



This normalization condition yields the following relation among the coefficients A, B of Eq. (fT7|) : 

\A\ 2 -\B\ 2 = -(-kr 1 )^a^. 



(20) 



(21) 



The solution associated with the adiabatic vacuum for k — ► oo corresponds to choosing A — {—ki]) 2 (p- 1 ) a 1 / 2 (— Zip) 1 / 2 
and B — 0. Thus one obtains the following adiabatic solution: 



which, in the proper time, becomes 
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ADIABATIC SUBTRACTION 



We now compute the integrals obtained on taking the vacuum expectation values of the relevant operators. Let us 
consider the correlator (ip 2 ) (we will link this to the energy-momentum tensor in the next section) 



<^ 2 > = 



1 



k\>l 



dk\ipk\ 



(24) 



where I = CHi is an infrared cut-off related to the begin ning of inflation [24], [25| . Using solution (j23|) we obtain a UV 
divergence, so we employ, as in the previous papers [lj, [l8|], dimensional regularization to treat the UV behaviour, 
therefore the integrands will be in 3 dimensions and the integration measure analytically continued in d dimensions. 
Subsequently an adiabatic subtraction is performed in order to obtain the renormalizcd quantities. The (ip 2 ) becomes 
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On using (B.3) of Appendix B of [18| with a. — d — 1 one obtains: 
1 



(^) 



d _ 3 r(3/2) r(-i/2)r(-i/2 + v ) 



2H; 

T{d/2 + v 




T{d/2) T(3/2 + ^) r(l - d/2)T{\ - d/2 + v) 
+0(d-3) 

where to avoid singularities in the analytic continuation we have to consider p > 3, this is not a problem because, for 
power-law inflation, the range p < 60 is disfavored at 2 a (see, for example, [261]). 
The adiabatic fourth order is (see Appendix and [HI]): 
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and the resulting renormalized quantity is 
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So the leading behaviour of the renormalized correlator (which is essentially given by the infrared contribution of 
((/? 2 ), namely the first term of Eq. (|26p ). during the inflationary era, with p >> 1, is: 



(^ 2 )ren 



V tt2 



(29) 



the same result is obtained in [13] for a massless minimally coupled test scalar field in a power-law model of inflation. 
In fact this test scalar field satisfies the same equation of motion of the field fluctuation tp (see Eq.(fT5j)'). 

Let us comment on the possible gauge dependence of this result. Starting from the general metric ([4]) and considering 
the first order gauge transformation (see, for example, (28[) the gauge invariant Mukhanov variable Q is given by 

ip+ jjijji and in the UCG this gauge invariant variable, as already stated, coincides with the scalar field fluctuation ip. 
Thus, in this particular gauge, ip 2 coincides with Q 2 which is a gauge invariant quantity. If we choose a different gauge 
such as, for example, the longitudinal one, the Mukhanov variable is no longer equal to the scalar field fluctuation tp 
and, while (Q 2 ) takes the same value as in the UCG, (tp 2 ) takes a different value. More generally for any quantity, 
in a given gauge, one can always choose a gauge invariant quantity which coincides with it in that particular gauge. 
The difficulty is that, except for some particular cases (such as ip 2 in the UCG), it may be difficult to understand the 
physical meaning of a general gauge invariant quantity. This appeares to be the case for the EMT, which is considered 
in the next section. A general gauge invariant variable connected with it in the UCG may not coincide with the value 
of the EMT in another gauge. So for the EMT the problem of the gauge dependence of its value does not appear to 
be solved. 



VI. APPROACHES TO THE BACK-REACTION 



One may follow different approaches in order to study the backreaction effects due to cosmological fluctuations. 
The main two methods in order to tackle this issue can be described in the following way. 
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The first consists of considering only first order perturbations, then imposing, to first order, the energy and mo- 
mentum constraints and finally defining an effective EMT by including all the quadratic terms present in the Einstein 
equations. Subsequently one averages everything over the quantum vacuum and first order quantities disappear. One 
finds that the effective EMT which appears in the averaged Einstein equations is modified by the back-reaction. 

The second approach is related to the standard perturbation analysis of the Einstein equations up to second order. 
In this case we impose the energy and momentum constraints and study the inflaton equation of motion perturbatively 
up to second order. One does not define any modified EMT but directly studies any observable averaged over the 
quantum vacuum in this framework. 



A. The Energy-Momentum Tensor of Cosmological Fluctuations 

On following the first approach the effective EMT of cosmological fluctuations is given by 

^fl _ /^biquadratic ][/[ 2 ^f/i quadratic (30) 



This method of considering the EMT of gravitational fluctuations is treated in textbooks [29| and has also been 

previously used in 0, H E3, ISSl . In this scheme one considers the modified Einstein equations M^Gp ^ = T„ + (rjf) 
which therefore includes back-reaction effects. 

For a generic potential, with an effective mass different from zero, the leading terms in the slow-roll parameter 
of the energy density are (see 



2 <^>(l-^). 



Analogously the average pressure is: 



(rj) = P 5) ~ 5} ( --^V) - 6H( V 2 ) )~~e8) (32) 



We now restrict the analysis to a power-law potential where r\ = 2e. 

The leading behaviour of the renormalized correlator is given by the infrared part of (ip 2 ) (the adiabatic part does 
not contribute to this leading value) so on using Eq. (f3~ljl we can obtain, by Eq. (|29|) . the leading value of the energy, 
without repeating the adiabatic subtraction, as 

eren ~ ~8^2^ 2 ^ ' fity 

Such a value leads to a negative contribution in the right hand side of the Einstein equation of the order of 
H 2 H 2 /Mp,, so we have that the average expansion rate appears to be decreased by the back-reaction of cosmological 
fluctuations. The importance of back-reaction is therefore related to the ratio Hf /M 2 l7 if inflation starts at a Planckian 
energy then back-reaction during inflation cannot be neglected as expected. We feel, however, that it is important 
to underline that for energy near to the Planck value the semi-classical approach to the quantum gravity may be 
questionable. 

We want to stress that, in this UCG, the back-reaction effects for a model of power-law inflation are very differents 
from the effects for a chaotic model ^(j) 2 found in [l3| and given by a negative effective energy of the order of 
Hf /H 2 . While in the case under consideration the effective energy has the same time behaviour as the background 
energy 3M^H 2 , in the chaotic case it grows with respect to the background value and a negligible back-reaction at 
the beginning of inflation can become important before the end of inflation (see [HI]). This is in agreament with 
what Abramo and Woodard state in [ll| on using the longitudinal gauge for the calculation of the backreaction in 
a power-law model of inflation, namely that the shape of the inflation potential can have an enormous impact on 
back-reaction. 

For the sake of completenes we have to repeat that the problem of the gauge dependence of this effective energy as 
yet does not appear to be solved. 
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B. Back- Reaction on the Geometry 



In the perturbative approach to the Einstein equations any back-reaction effect is analyzed by evaluating pertur- 
batively quantities which characterize the geometry (as, for example, the expansion scalar 0). As stated in Sec. [V] 
the leading effect in the renormalized quantities, for models with an effective mass different from zero, is given by the 
infrared part of the non renormalized integrals, so we shall use the long-wavelength approximation in the treatment 
of the second order scalar perturbations (see [H, HH). In the long- wavelength limit, from the first order equation 

of motion Eqs.© and ([10)) , one obtains ip ~ Cjjj where C is a constant in time, thus in this limit the non-local 
spatial contribution s (Eq. (|13|) ) has an ordinary behavior given, for the isotropic case, by cptp / (AM^H) . In general on 
using this limit on the right hand side of Eq. (fT2)) one obtains the following result (valid to all orders in the slow-roll 
parameters): 




(34) 



For the case under consideration, namely power-law inflation, 4s is a constant and, consequently, the inhomogeneous 
term of Eq. (|12[) is equal to zero in the long- wavelength limit. A coupling between the modes of different orders in 
perturbation theory could then occur only in the ultraviolet. Thus, in the long-wavelength approximation, the equation 
of motion for the second order field fluctuation is the same as that for the first order one, and one can write in 
the infrared an equation for a total fluctuation (p = ip + p 1 - 2 ^ which is to be constrained by the initial condition. This 
is equivalent to saying that since the first and the second order scalar perturbations ip and p>^ are decoupled, the 
second order contribution does not give us any new information, with respect to the first order, and it can be seen as 
a renormalization of the first order one. This is the only useful information which can be gained by the perturbative 
approach. We feel that the significance of the above considerations and results deserves further investigations. 



VII. DISCUSSION AND CONCLUSIONS 



The renormalized EMT of cosmological fluctuations for a power-law model of inflation has been studied. As in our 
previous work [13| we have self-consistently taken into account the gravitational fluctuations accompanying the scalar 
field fluctuations in a UCG. 

We find that the renormalized EMT of cosmological fluctuations during slow-rollover carries negative energy density 
(due to the inclusion of gravitational fluctuations [l3[) with a de Sitter like equation of state to leading order. The 
average expansion rate appears to be decreased by the back-reaction of cosmological fluctuations and the leading value 
of the effective energy is of the order of H 2 H 2 , so the back-reaction effects cannot be neglected if inflation started at 
Planckian energies, i.e., at Hi ~ M p \. 

We stress that the back-reaction effects for this model arc very diffcrcnts from the effects for a chaotic model 
([HI) gi ven by a negative effective energy of the order of Hf / H 2 . While in the case under consideration the effective 
energy has the same time behaviour as the background energy SM 2 X H 2 , in the chaotic case it grows with respect to the 
background value and also a negligible back-reaction at the beginning of inflation can become important before the 
end of inflation (see [HI). The back-reaction effects in inflation are strongly dependent on the shape of the potential 
(see [ll| for similar consideration in longitudinal gauge). 

In this model of power-law inflation and in the long- wavelength approximation we also find that the first and second 
order scalar perturbations <p and <p>^ are decoupled and they have the same equation of motion, the second order can 
then be seen just as a renormalization of the first order with no other physical information. 

To conclude we find in the Appendix the fourth order adiabatic expansion for a general potential V(<p). Let us note 
that to obtain this general fourth order adiabatic expansion a crucial point is to distinguish between the different 
adiabatic orders, the effective mass Va,a obtained by a double field derivative is of adiabatic order while the quantities 
obtained by deriving n-times with respect to the time are of adiabatic order n, independently of their value. 
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VIII. APPENDIX: THE FOURTH ORDER ADIABATIC EXPANSION 



In order to remove the divergences which appear in the integrated quantities as poles in the T functions, we shall 
employ the method of adiabatic subtraction [3ll [32J. Such a method consists in replacing our functions with an 
expansion in powers of derivatives of the logarithm of the scale factor. This expansion coincides with the adiabatic 
expansion introduced by Lewis in [33[ for a time dependent oscillator. 

Usually it is more convenient to formulate the adiabatic expansion by using the modulus of mode functions Xk — 
l^k/ \/2| and the conformal time r\ [32| [dr\ = dt/a). We follow this procedure and write an expansion in derivatives 
with respect to the conformal time (denoted by ') for Xk- We then go back to cosmic time and insert the expansion 
in the expectation values we wish to compute. Adiabatic expansions in cosmic time and conformal time lead to 
equivalent results, because of the explicit covariance under time reparametrization (34[. 

The variable Xk satisfies the following Pinney equation: 



Xk + 3Hik 



k 2 H 
1? V « + 2 ( 3H+ H 



Xk 



a 6 xl ' 



(35) 



Following [13l . [18| we rewrite Eq. (|35p in conformal time in the following way: 

[ax k )" + Q 2 k (ax k ) ' 



(ax k ) s 



(36) 



where 



Ql = k 2 + o 2 U 



-a 2 R 



(37) 



and R is: 



/ a' 2 a" 2 
R = R-6\ -4— -2-=-^ 



(38) 



with R = 6^3- the Ricci curvature. The key point for this calculation is that is of adiabatic order 0, while R is 
of adiabatic order 2, independently of their value. From Eqs. (|36I37I38|) one obtains the expansion for Xk up to the 
fourth adiabatic order: 



„( 4 ) 



1 1 



i 



' - 1 1 ~ I 62 
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where Qk is defined in Eq. (|37[) and £2 , £4 are given by: 



€2 ~ ~2W + aW 



m k . 1 1 » 

" 4 " 4«f 2 " 
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The solution in Eq. (|3"!)|) must be expanded again since R is of adiabatic order 2. Therefore x^\ for a general 
potential, is: 
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1/2 v l/2\ 
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where c = a 2 and 
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